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ABSTRACT 



This is a continuation of the work of assembling and testing 
the equipment necessary to conduct gairana^gamma angular correlation 
experimeotSj started by Fort A. Verser, Jr* at the JS Naval Post- 
graduate School (1). 

No substantial changes were made in the equipment^ bat rather 
a rerun of the Cobalt-60 experiment was made to confirm the results 
obtained in the previous work. A comparison of the previous work 
confirms that there is a definite distortion in the shape of the 
observed correlation function. 

The statistical analysis of the data and the solid angle 
corrections required for evaluation of the results have been pro- 
grammed for calculation by use of the 1604 Digital Computer* This 
paper outlines these programs in detail. 

The author wishes to express his appreciation for the assist- 
ance given him by Professor Harry E. Handler and tor his guidance 
throughout this investigation. Appreciation is expressed to 
Professor Edmund A. Milne for his advice and commentSo Captain T. 
R. Abernathy^ USMC, was very helpful in providing valuable assist- 
ance and discussions on various aspects of the computer programing. 
The cooperation and assistance of Mr. Mervyn C* Brillhart in main- 
taining and testing the electronic equipment was indispensable* 
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CHAPTER I 
INTRODUCTION 



The decay of Cobalt-'60 to an excited state of Nlckel“60 and 
the angular correlation of the resulting gaimna- gamma cascade have 
been studied extensively experimentally (25638595 and others) and 
the results of the correlation measurements agree very well with 
the theoretical correlation for the decay scheme shown in Figure 1. 
Consequently 3 the decay of Cobalt«60 is now used as one standard 
for the testing of correlation equipment. 

Angular correlation equipment was assembled at the US Naval 
Postgraduate School by Verser in 1960 for the investigation of the 
properties of low«lying excited nuclear states. His testing of the 
operating characteristics of the equipment by the angular corre^ 
lation of the cascading gamma rays resulting from the decay of 
Cobalt 60 yielded an observed correlation curve whose shape was not 
in agreement with the theoretical curve. This work was Initiated 
in order to attempt to verify Verser *s results. 

In order to facilitate the statistical analysis of the data^ 
a program has been worked out for the use of the 1604 Digital 
Computer, The program follows the procedures outlined by MoE, 

Rose (2) and is contained completely in Appendix 1, 

An additional program has been written to compute the solid 
angle correctionSs which must be applied to the theoretical cuive^ 
by the method of M,E, Rose (2), This progam is outlined in Ap- 
pendix II, 

The results of this work agreed within statistical expectation 
with those obtained by Verser 5 indicating that the imexpected shape 
of the curve is apparently real and not an instrumental effect. 
These results are given in Chapter VI, 
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CHAPTER II 



THEORETICAL CONSIDERATIONS 

The deterini nation of the angular momentuni quantum numbers^ 
or spins 5 of low-lying short^^lived excited nuclear states by 
the nuclear spectroscopic tool, known as angular correlation of 
successively emitted radiations, is quite common practice today. 

Only within recent years, because of the vast advancement of 
electronics, has practical utilization of this tool become 
possible • 

Since nuclei, under ordinary conditions, are randomly ori» 
ented in space, it is not possible to observe a radiation pattern. 
This is due to the fact that the probability of emission of a 
radiation by an excited nucleus, depends on the angle between the 
nuclear spin axis and the direction of the emission. Therefore, 
in order to observe a radiation pattern, it is first necessary to 
orient the nuclei. 

One method of orienting certain nuclei consists of placing the 
sample in an electric field gradient or a strong magnetic field at 
a very low temperature. Another method, the one used in the present 
consideration, consists of selecting nuclei whose nuclear states 
have spins lying in preferred directions. This happens to be the 
case for the intermediate state if a nucleus decays by successive 
emission of two radiations. 

In an angular correlation experiment the first radiation is 
observed in a fixed direction. This establishes an axis to which 
the direction of emission of the second radiation, originating from 
the state formed by the first, can be referred. The second radiation 
has preferred angles of emission with respect to the direction of 
the first. 

In order that the correlation exhibit maximum anisotropy, the 
angular momentum vector of the intermediate state must not change 
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direction significantly before the second radiation occurs » The mean 

“8 

life of the intermediate state must be less than about 10 seconds 5 

the typical precession period of the nuclear spin about local per- 
turbing fields. The angular correlation function is determined by 
the nature of the radiations and by the spins of the states involved, 
and thus its measurement may lead to an unambiguous choice of spins 
for the states. 

The theoretical expressions for the correlation functions have 
been worked out for a large number of cases of interest (7) . For 
the particular case of a gamma-gamma cascade, the correlation 
function W(0) is 



n n 

max^ max^ 

W(0) = / , A P (Cos 0) = / , B Cos"e (1) 

n = 0 n = 0 

where P is the Legendre Polynomial of even order n, n is the 
n max 

smallest even integer of the set of three numbers consisting of 
twice the multipole order of each of the two gamma rays 5 and 0 is 
the angle between the directions of emission of the two gamma rays. 

The theoretical correlation function for the cascade in Nickel-bOj, 
assuming point detectors and a point source, is 

W(0) = 1 + 0.1020 P^CCos 0) + 0.0091 P^(Cos 0) (2) 

with an anisotropy 

R ^ W(180°) - W_(90°2_ 0.1667 (3) 

W(90°) 



Many workers, after correcting these theoretical results for the 
effect of non-point detectors have found good agreement between them 
and their experimental correlation results. 
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CHAPTER III 

DESCRIPTION OF EQUIPMENT 



The equipment needed for the gamma-gamma correlation 
measurements must furnish means for detecting the radiation^^ 
pulse-height amplification with selection or discrimination, 
coincidence analysis, and adequate recording of the counts in 
the various counting channels. A mechanical system must also be 
provided for mounting the detectors and changing their relative 
positions and for holding and positioning the radioactive source. 

The apparatus used for this work was that assembled by Verser 
and described in detail in his report. Therefore, only a brief 
description of the equipment and of the minor alterations per- 
formed on it will be presented. 

Figure 2 is a photograph of the mechanical setup. The origi- 
nal base of the correlation table has been replaced, the new base 
is of solid brass with a bearing mounted shaft to enable accurate 
manual detector positioning. 

The electronic equipment (Figure 3) has the functions of 
detection of the radiation, pulse-height amplification and dis- 
crimination, coincidence analysis, and recording the number of 
counts. The only change in the electronic equipment was in the 
pulse-shaping networks used in the initial calibration of the 
apparatus. The pulse shapers were rebuilt to give a better simu- 
lation between the test pulses and the actual pulses encountered 
under operating conditions. 



6 









I 






IM 

Z 

z 

< 

X 

o 




q; < 

o 

H a: 
o < 

LU LlI 

»- 2 
UJ — 
Q -I 



(T 

O 

5 



IT 



a: ^ 



q! < 



(T 

UJ 



.T Q -I 



CL 

Q. 

3 

CO 

(T 

UJ 

o 

0 . 



3 CO 



< UJ 

o o 



H 

2 

UJ 

cr 

UJ 

u. 

o 



UJ 

CD 

3 



< 

!3 

o 

> 



go 

e> I 



I 



I I 

Ul S <2 ^ CO < 
Q < Q CO 0- O 



>“ 

UJ 



• s 



SCHEMATIC DIAGRAM OF THE ELECTRONIC SYSTEM 



CHAPTER IV 

EXPERIMENTAL ANALYSIS 



The angular correlation function W(0) (Equation 1) nmst be 
expressed in terms of experimentally-observable quantities. The 
equipment records the total number of coincidences as a function of 
the angle between the detectors, and the total number of radiations 
captured in each channel. C^, the total number of coincidences in 
a counting period, is given by 



C — C 4 - C + C- 
t g a 



(^) 



where 



C = the number of genuine coincidences due to the two 
^ cascading radiations, 

C - the number of accidental coincidences due to the 

Si 

finite resolving time of the electronic equipment, 
and 

= the number of background coincidences. 



The half“life of the source is long compared to the time required 
to collect the necessary data; therefore, equation (4) can be re- 
arranged and rewritten in terms of counting rates. 



C 

g 




“ C 

a 




(5) 



C can also be 
g 



expressed in terms of other parameters: 



C 

g 



(C a P oc P + £ a 9 £ a P n q ^ w^Q ) (6) 



12 12 12 21 21 ' 21 ^ 1 2^12 0 



where 

N = the absolute source strength 
o 

0 = the fraction of nuclei decaying by the emission of 

^ gam.ma ray ’Vj" and gamma ray ’’k" in a cascade 
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“^i = the fractional solid angle subtended by the crystal 
in Channel "i" 

= the fraction of garnina rays *'j*' which leave the source 
in the right direction to enter crystal "i" that are 
not absorbed before reaching the crystal 
£ = the total efficiency of the crystal in Channel "i" 
for gamma ray 

= the fraction of the number of gamma rays captured 

by the crystal in Channel *'i" whose pulse heights are 
accepted by the energy discrimination of Channel " i” 
W(0 ) = the angular correlation function at an angle 0 
averaged over the solid angles used, 

the counting rate in Channel *'i", can also be expressed as a 
function of certain of these same parameters. 



N. = A, 0.N 

1 11 




.a. . e . . 

IJ 1 J ^ IJ 



+ N 



ib 



(7) 



where 

0^ = the number of gamma rays "i'* per disintegration 
Nib ~ background counting rate in channel *'i” 



From Equations (5) and (6) we get for the correlation function^ 




where f is the appropriate function of the £*s, 0 :''s 5 and 0*s. 
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By combining Equations (7) and (8), we get for the correlation 
f unc t i on 



W(0 ) = 




(9) 



where g is a function of the s 0’s. 



In this form W(0 ) is relatively insensitive to equipment 
m 

Therefore;, it is only necessary to calculate the modified 
function given by 



^ W(0 ) 

_ m 

W(0 ) = 

m N 

g o 



r. 



(C -C -a ) 
^ t a 






[^1 



J 



m 



stabil Ity • 
correlation 



( 10 ) 



The correlation runs were done with the movable detector at 
seven angles; 90°, 105°, 120°, 135°, 150°, 165°, and 180°. A 
minimum of six two-hour counting periods were completed at each 
angle, and the angular order in which they were taken was varied 
to minimize any systematic f luctuations o Each run yielded single 
channel counting rates of approximately 1200 counts per second and 
a total coincidence rate of approximately Id counts per second. 

Before the experimental results can be compared to theory, the 
data must be checked and corrected, if necessary, for the various 
deviations from the ideal arrangement, i.e,, centered point source, 
no spurious coincidences due to the coincidence resolving time, no 
scattering, no background, and perfect stability of the electronic 
equipment. In addition, the theoretical correlation function must 
be corrected for non-zero solid angle detectors o The methods for 
accounting for these effects are discussed in the following para- 
gr aphs o 



11 



J 



i 



A. SOURCE SIZE AND POSITIONING 



As a test of the positioning of the source in the center of the. 
correlation table, the integral counting rate of a portion of the 
Cobalt“60 spectrum was measured in each channel as a function of the 
angle between the counters. This was done by Verser and checked again 
after the base of the correlation table had been changed. Within sta- 
tistical expectations, the counting rates were constant. The source has 
a volume of about 0.001 cubic inches and can be regarded as a point source 
at the distance at which the detectors were placed from it. 

B . SCATTERING 



Compton scattering occurring outside the detectors can give rise to 
unwanted coincidences. These in general will tend to smear out the mea- 
sured correlation function. This problem is minimized by use of scintil- 
lation crystals as detectors and accepting only the full-energy peaks of 
the desired gamma rays. The method for determining the pulse height se- 
lection limits has been discussed in detail by Verser, 

C. BACKGROUND CORRECTION 



With the discriminators at the same settings as were used during the 
correlation runs, background runs were made before and after the corre- 
lation runs. The background coincidence counting rate was negligibly small 
and the background rate of each channel was such as to yield total back- 
ground counts of the order of 0.4 percent of the total counts in each 
channel during the correlation runs. These background rates were used as 
corrections to the gross counting rates obtained during correlation runs. 

D. DETERMINATION OF GENUINE COINCIDENCE RATE 

Since the background coincidence counting rate proved experimentally 
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j 



to be negligible. Equation (5) becomes 



C = C “ C 
g t a 



For the total counting time, T, the total number of coincidences is C^T 
and the standard deviation in the. total coincidence rate is 




Therefore, the standard deviation of the. genuine coincidence rate is 



2 

(T 

g 






( 12 ) 



The accidental coincidence rate is given by 



C 

a 



2 ^^ 
r 1 





o 



(13) 



where. 



- Coincidence resolving time of the coincidence analyzer. 

In order to minimize, the accidental coincidence rate, the. source 
strength must be kept as small as possible. However, if the source 
strength is decreased, the counting time must be increased to obtain a 
specified precision in the counting rate. Thus, a compromise must be 
reached between precision and counting time. 

The smallest resolving time obtainable, with the present equipment is 
experimentally determined to be 

= (0,1820 + 0.0016) microseconds; 

this res il ted in an accidental coincidence, rate that was about 307o of the 
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total coincidence rate. The approximate source strength of the water 
solution of CoCl^ used in this experiment was 0.05 millicuries of Cobalt‘-50, 



E. REFERENCE COUNTING RATE 



With the discriminator in each channel adjusted to accept only the 
full^energy peaks of the desired gamma rays, the integral counting rate 
of the Cobalt“60 spectrum in each channel was determined o This gave a 
reference counting rate for each channel. During the course of the corre^ 
lation runs, either the discriminator settings or the gains of the linear 
amplifiers were adjusted before each run so that, the integral counting 
rate agreed within + 5 cps with these reference rates (Table 1). 

TABLE 1 

REFERENCE COUNTING RATES 

N^ = (1085.26 + 0.57) cps 

N^ - (1117.32 + 0.59) cps 

As a measure of the gain fluctuations during each counting run, the 
integral counting rate of each channel was measured after the run using 
the same discriminator setting that was employed during the run. The 
difference between this counting rate and the reference counting for that 
channel was determined. If the variation exceeded + 20 cps in either 
channel, the run was discarded. 

This procedure resulted in about 307o of the runs being discarded 
because they fell outside of the maximum acceptable fluctuation. This 
indicates that the electronic system is not as stable as it should be. A 
check of the line voltage revealed substantial fluctuations during any 2^« 
hour period. The temperature of the room was not controlled in any manner 
and varied considerably. These could have been the major causes of the 
large counting rate fluctuations since the discriminators are highly sensi- 
tive to both temperature and voltage variations. 
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F. SOLID ANGLE CORRECTIONS 



The solid angle correction which must be applied to the theoretical 
correlation function to account for the non-zero solid angle of the de- 
tectors and enable a comparison to be made with the measured function is 
described by Rose (2) . The correction for two cylindrical crystals whose 
axes intersect at the source involves the numerical evaluation of integrals 
of the form 



n, 1 



_r X (0) 

(Cos 0) ^ 1 -e J Sin 0d0 



where 



(14) 



- the full-energy absorption coefficient of the detector in 
Channel "i” for the gamma ray; 

^ = the half-angle subtended by the front face of the crystal 
in Channel "i”; 

-1 -r . 

X.(0) = t.Sec0 for 0^0:^ Tan i__ ; 

^ ^ h . + t . 

1 1 



X^(0) = r Xsc0-h^Sec0 



for Tan 



-1 



J 50 

+ t . 1 



where 

h. = the distance from the source to the crystal in Channel "i": 
1 

t. ^ the thickness of the crystal in Channel ”i”^ 

1 

r. = the radius of the crystal in Channel "i"o 
1 



The attenuation factors, Q^, are functions of these integrals: 



’n.l 




I „ 
n, 2 


'o.l 




I 0 
0,2 






The corrected correlation function becomes 



( 15 ) 



w(e) = 



I 



max 
1 



Q A P (CosQ) 
n n n 



n=o 
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CHAPTER V 

EXPERIMENTAL RESULTS 



t It. 

The modified correlation function for the k run at anele 0* aiven 



by Equation (10) is 



K («») =■ 



- c - c, 

tab 



A 



(Nj-Nib) (K^-Njj,) 



This means that the seven counting rates C , C . C . N, , N,, , and 

^ t’ a" b" 1" Ib"^ 2’ 2b 

must be determined for each run at each angle and the standard deviation 
of each rate calculated. 

The background counting rates and for each run were expert^ 

mentally determined. The single-channel counting rates corrected for 

scaling losses yielded and The accidental counting rates were cal 

culated for each run from Equation (13). The total coincidence counting 

rate, C^, was experimentally determined. The background coincidence rate 

was negligible throughout the experiment. 

The experimental value of the modified correlation function W/ (9*) 

Km 

and its standard deviation were calculated for each run at each angle. 

From these the weighted mean W* (9*) and its standard deviation were com-^ 

m 

puted for each angle. 

The weighted mean of the measurements taken at each angle was deter- 
mined by using the inverse square of the standard deviation of each mea- 
surement as a weighting factor: 




( 18 ) 



where 



N/ = number of runs at the angle 0^ . 
m ^ m 
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The standard deviation of the weighted mean 
the individual standard deviations: 




/ 



(W") was calculated from 



(19) 



A comparison between the standard deviations of the means calculated 
from Equation (19) and those calculated from 



cr^(w') 




( 1 


) -W(0 


1 » 


N.( 


N„-/) 





( 20 ) 



revealed no significant differences. Therefore the standard deviation 
obtained from Equation (19) were used in subsequent calculations and the 
method discussed in paragraph E Chapter IV for treating gain instabilities 
is realistic. 

The values of the weighted means and their standard deviations are 
tabulated in Table 2. 



TABLE 2 



The Weighted Means W* (O^) 



0. (degrees) 
m 


Weighted 

Mean 


Standard 

Deviation 


90 


5046 


45 


105 


4925 


33 


120 


5109 


21 


135 


5211 


36 


150 


5556 


35 


1 165 


5667 


24 


1 180 


5907 


43 
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A curve of the form W’(9) ^j^(Cos0) was fitted to the exper i 

mental by the method of least squares as outlined by Rose (2). 

This analysis was programmed to be carried out by the 1604 computer. (See 
Appendix II for details). 

The values of A* are given by the matrix equation 



A* 



T 

A \ 



1 . 



•I Mil 



( 21 ) 



where 



U) 



I|a| 



T 

A \)A 



= the 7x3 matrix consisting of the Legendre Polynomialsj 
P^, P^j and P^, evaluated at the seven angles (Figure 5A) o 

= the transposed A matrix 

= a diagonal matrix with elements consisting of the inverse 
squares of the standard deviations of the experimental 
W*(9^) (Figure 5B) 

= a one column matrix whose elements are the experimental 
values of W' (0^) (Figure 5C) 

= the one column matrix whose elements are respectively 
the least“square coefficients A^j A^ and A^ (Figure 5D) 

the inverse of the matrix obtained by myltiplying the 
three designated matrices together. The diagonal terms 
of this matrix are the squares of the standard deviations 
of the least-square coef f icients^ A** . 



The experimental values for the coefficients A® and their standard 

n 

deviations, normalized so that A* = 1. are listed in Table 3, The corre- 

o 

spending values from Verser's work are included for comparison. 
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TABLE 3 

Least Square Coefficients 



Ver ser 






This Work 


A' 1.0000 


4 - 


0.0034 


1.0000 + 0.0026 


o 








Al 0.0961 


± 


0.0055 


0.0971 0.0057 


2 








a! 0.0339 


± 


0.0071 


0.0231 4 0.0066 


A 
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MATRIX CONFIGURATIONS 



P^(Cos90°) 


P2(Cos90°) 


P^(Cos90°) 


P (Cosl05°) 
o 


P^CCosIOS®) 


P^(CosI05°) 


P^(Cosl20°) 


P2(Cos120°) 


P^(Cosl20°) 


P (Cosl35°) 
o 


P^CCosUS®) 


P^(Cosl35°) 


P^(Cosl50°) 


P^CCoslSO®) 


P^(CosI50°) 


P^(Cosl65°) 


P2(Cos165°) 


P^(Cosl65°) 


P^(Cosl80°) 


P^CCoslSO®) 


P^(Cosl80°) 



(A) 



CT"^(90°) 


0 


0 




0 


0 


0 


0 


0 


C/'^(105°) 


0 




0 


0 


0 


0 


0 


0 


(T“ 


^(120° 


) 0 


0 


0 


0 


0 


0 


0 




CT-' 


^(135°) 0 


0 


0 


0 


0 


0 




0 


(^"^(150°) 


0 


0 


0 


0 


0 




0 


0 


CT“^(165°) 


0 


0 


0 


0 




0 


0 


0 


cr-^iso” 








(B) 








W' (90°) 
W' (105°) 
W' (120°) 
W' (135°) 








A 

O 








W' (150°) 
W' (165°) 
W' (180°) 








(d: 


) 






(C) 

















F IGURE 5 
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A> Experimental Anisotropies 



The anisotropy is defined: 

R ^ W(180°) - W(90° ) (22) 

W(90°) 

The anisotropy for the least-square curve was calculated from the experi- 
mental correlation function. Since the anisotropy can be determined with 
much greater precision than can the least-square coefficients themselves^ 
it provides a more sensitive means for comparing the experimental and theo^ 
retical functions. The experimental anisotropy is 

R - 0.167 + 0.013 

The standard deviation was calculated by the formula presented by Klema 
and McGowan (6). The computer program for the calculation of R and its 
standard deviation is outlined in Appendix I. 
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CHAPTER VI 

RESULTS AND CONCLUSIONS 
A. COMPARISON BETWEEN EXPERIMENT AND THEORY 



The experimental and corrected theoretical functions for this experi- 
ment are: 

Experimental • 

W*(0) = (1.0000 i 0.0026) + (0.0971 + 0.0057) P^ 

4- (0.0231 + 0.0066) P^ 

Theoretica 1 : 

W(0) - 1.0000 + 0.0996 + 0.00838 P^ 

The experimental anisotropy was calculated to be R - 0.167 + 0.013 and 
is compared with Verser*s value of R = 0.171 j_ 0.027 and with the theo-- 
retical anisotropy corrected for the solid angle^ R = 0ol61, 

Figure 6 shows the experimental results and their standard deviations for 
this work and that of Verser. Also shown are the least-square correlation 
curves for both sets of data and the corrected theoretical correlation 
curve. All the curves have been normalized to unity at 0 = 90°. 

B. DISCUSSION AND CONCLUSIONS 

From Figure 6^ it can be seen that the curve obtained In this work 
agrees with that obtained by Verser. Both of these curves show an ap- 
parently real dip in the region from 105° to 135°. It is only a remote 
possibility that this dip can be attributed to statistical f luctaa tions^, 
and it is unlikely to be due to instrumental effects. 

Two conjectures of possible mechanisms are: 

(a) The source may have crystallized resulting in preferential 
alignment of the nuclei within the crystal due to Internal 
electric fields. 
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(b) The ferromagnetic properties of Cobalt could have led to 
a preferential alignment of the nuclei in the earth’s cnag^ 
netic field. 

The following series of experiments are suggested lor future work. 

(a) The results for measurements conducted for angles between 
180° and 270° should be compared with those already ob“ 
tained . 

(b) In order to check on the possible effects of orientation 
due to crystalline electric fields, the experiment should 
be repeated with the present source at various orientations 
relative to the fixed counter. 

(c) In order to check on the possibility of orientation by vnag- 
netic effects the source should be placed in a magnetic field 
of known direction and the measurements repeated, 

(d) A correlation curve using a powdered Cobait-60 source in- 
place of the present source should be obtained for comparison , 

These experiments would determine conclusively if any effects causing pre- 
ferential nuclear orientation are present. 

The following instrumentation improvements are suggested in order to 
improve the stability of the equipment. 

(a) A very stable power supply should be purchased to aid in 
preventing the line voltage fluctuations from interfering 
with the electronic system. 

(b) The effect of temperature variations on the electronic system 
can be minimized by placing the equipment in a temperature 
controlled room for conducting future work, 

(c) A slow coincidence unit whose resolving time is independent 
of the counting rate should be obtained. 

The precision of the results should be improved by the addition oi a 
fast coincidence channel working in conjunction with the slow coincidence 
unit to reduce the number of accidental coincidences. 
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Since the experimental curves using Cobalt-60 have not agreed with 
the theoretical curve, the experiments on Cobalt-60 should be continued 
until the source of the distrotion is ascertained » Perhaps the effect 
can be related to the physical properties of the CoCl^ source* 
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Least-Square Fit of Experimental Results 

The comparison of the experimental correlation function to the 
functional form of the theoretical correlation function given by Equation 

(1) is accomplished by the method of least-squares, as outlined by Rose 

(2) . The analysis requires the solution of Equation (20) for the values 

of A* . 
n 

The computer solution of the matrix equation, like the solid angle 
calculation, is carried out in floating point format « In addition to the 
solution for the values of the program also calculates the anisotropy 

by Equation (21), and its standard deviation by the formula presented by 
Klema and McGowan (6). A flow sheet outline of the program is given in 
Figure (7) o In order to provide a certain amount of flexibility to the 
program, the angles used have been left as a parameter to be supplied. 

The program as written is, however, limited to seven different angular 
positions. The parameters needed for the program may be placed anywhere 
in the computer and their addresses supplied to the program by placing 
them in the B registers in the following order: 

B1 = weighted means of the correlation function in ascending angular 
order , 

B2 the standard deveations of the weighted means in the same order 
as above. 

B3 = the values of the angles used in ascending order. 

The program is available for future use on both bioctal tape in 
machine language, and IBM cards in assembly language. The bioctal tape 
contains the following: 

1. The main program for calculation of A^ 

2. The program of normalization of the A^ 

3. The program for calculation of the resulting correlation function 

4. The normalization of the correlation function 

5. The calculation of the anisotropy 
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ANISOTROPY magnetic TAPE 
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6« The calculation of the standard deviation for the anisotropy 
7, The subroutines 

a. Floating Point Sine-Cosine 

b. Floating Point Square Root 

Co Floating Point Decimal Output. 

The results of the computations are dumped on magnetic tape unit (4) 
in the following sequence: 

Lines 1 thru 3 The matrix inverse 

Lines 4 thru 6 Check answer; included in the program is a check 



of the inverse which consists of multiplying the 
inverse by the original matrix 



Line 


7 




The 


values of A* , A* and A* 
o Z 4 


Line 


8 




The 


standard deviations of the A* 

n 


Line 


9 thru 


13 


The 


resulting correlation functions 


Line 


14 and 


15 


The 


normalized correlation functions 


Line 


16 




The 


normalized A* 
n 


Line 


17 




The 


normalized standard deviations of the A* 

n 


Line 


18 




The 


anisotropy 


Line 


19 




The 


square of the standard deviation of the anisotropy 
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STATISTICAL ANALYSIS PROGRAM 
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APPENDIX II 



Solid Angle Correction 

The solid angle correction, which must be applied to the theoretical 
correlation function to enable a comparison to be made with the experi- 
mental correlation function, is described by Rose (1) o The correction 
involves the numerical evaluation of the integrals given by Equation (14): 

1^. - P^(Cos0) Sin0 d0 

o 

and the subsequent solutions of Equations (15) and (16). 

The computer solution for this correction was carried out in floating 
point format. The numerical integration was computed using an interval 
of 0.001 radians. A flow sheet outline of the program is given in Figure 
(8) o Ihe program must be supplied with the parameters h, r, and t 

same in both channels (see Table 4) thus enabling solid angle corrections 
to be made for any experimental set“up to which these corrections may be 
applicable. These parameters may be placed anywhere in the computer and 
their addresses supplied to the program by placing them in the B registers 
in the following order 

Bl = r (address of the crystal radius in centimeters) 

B2 t (address of the crystal thickness in centimeters) 

B3 - h (address of the crystal to source distance in centimeters) 

B4 - T (the value of the full-energy absorption coefficient of the 
detector for the gamma ray in centimeters ^) 

The values of the parameters used in the calculation of the corrected 
angular correlation coefficients for the theoretical correlation function 
and the results of the calculations are listed in Table 4. 

The program is available for future use on bi octal tape in machine 
language and on IBM cards in assembly language. The bioctal tape contains 
all the necessary subroutines to enable a complete computation, i.e. 

1. The main program for solution of Equation (14), 

2. Ihe subroutines 
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SOLID ANGLE CORRECTION FLOW CHART 



a. Floating Point Sine-Cosine 
bo Floating Point Exponential 
Co Floating Point Decimal Output o 

The results of the computations are dumped on magnetic tape unit (4) 
in the following sequence: 

Line 1 



Line 2 



I 



2 



Line 3 



1 



4 



Line 4 




Line 5 




Line 6 




Line 7 




Line 8 




Line 9 




6e 





TABLE 4 


SOLID ANGLE 


CORRECTION PARAMETERS 


h (cm) 


18.0 


t (cm) 


2.54 


r (cm) 

T (cm"^) 


2.54 

0.098 


^0 


0.0079 


^2 


0.0078 


'a 


0.0076 


^0 


0.9882 


■^4 

'o 


0.9609 


^2 


0.9765 


^4 


0.9234 




0.0996 


^4^4 


0.0084 
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SOLID ANGLE CORRECTION PROGRAM 





ORG 


12000 




REM 


BURTON J CCSTWAY APRIL I 96 I 




REM 


SOLID ANGLE CORRECTIONS 




REM 


OF THE THEORETICAL CORRELATICN 




REM 


CURVE. MUST BE SUPPLIED 




REM 


WITH PARAMETER ADDRESSES AS 




REM 


FOLLOWS. 




REM 


B1 ADDRESS OF XTAL RADIUS. 




REM 


B2 ADDRESS OF XTAL THICII'IESS. 




REM 


B3 ADDRESS OF CRYSTAL TO 




REM 


SOURCE DISTANCE. 




REM 


EL address of PARAIIETER TAU. 


HOLD 


EQU 


5500 


TiCMP 


ECiU 


6500 


WORK2 


EQU 


7500 


work4 


EQU 


10500 


HANS 


ECU 


11500 


SOIJI 


ECU 


12500 


JAYO 


ECU 


13500 


JAY2 


EQU 


15504 


JAY4 


EQU 


13510 


JAY20 


ECU 


13514 



64 



12000 

12001 

12002 

12003 

12004 

12005 

12006 

12007 

12010 

12011 









jay 4 o 


EQU 




13520 








EXPFL 


EQU 




60600 








SINFL 


EGjU 




60000 


12 


1 


00000 




IDA 


1 


0 


20 


0 


12163 




STA 


0 


RAD 


12 


2 


00000 




IDA 


2 


0 


20 


0 


12164 




SPA 


0 


THK 


12 


3 


00000 




IDA 


3 


0 


20 


0 


12165 




STA 


0 


HOT 


12 


h 


00000 




IDA 


4 


0 


20 


0 


12166 




STA 


0 


TAU 


12 


0 


12171 




IDA 


0 


ZRO 


20 


5 


13500 




STA 


5 


JMO 


54 


5 


00030 




ISK 


5 


30 


75 


0 


12004 




SU 


0 


/=! 


12 


0 


12171 




IDA 


0 


ZRO 


20 


0 


12151 




STA 


0 


DELTA 


50 


1 


00000 




ENI 


1 


0 


50 


2 


00000 




Eia 


2 


0 


50 


3 


00000 




EKl 


5 


0 


50 


4 


00000 




ENI 


4 


C 


50 


5 


00000 




ENI 


5 


0 


50 


6 


00000 




ENI 


6 


0 



65 





12 


0 


12163 


LDA 


KJ 


ILD 


^ORMATK'!: h/U 




35 


0 


12165 


FDV 


0 


HOT 




12013 


20 


0 


12172 


STA 


0 


PLEAS 






50 


0 


00000 


ENI 


c 


0 




12014 


12 


0 


12165 


LDA 


(; 


HOT 


FOWJSKM Ot 




50 


0 


12164 


jiiV-D 


0 


THK 




12ul5 


20 


0 


12175 


b'TA 


0 


a?CE 






12 


Cj 


12165 


ld; 


0 


K;j) 




12( >16 


55 


'J 


12175 


FDV 


0 


a-JCE 






20 


0 


12175 


STA 


0 


aiGE 




12017 


12 


0 


12151 Kiii 


U)P 


0 


DELTA 


rXfr -f.n 




50 


0 


00000 


LLI 


0 


0 


1 l.Vu'V til'' 


12u20 


75 


4 


60OCI 


SIJ 


4 


SIKE'L+1 


000124] 




50 


0 


000(J0 




0 


0 




111 i21 


50 


0 


ooooo 


au 


.) 


0 






50 


0 


coooo 


•aNi 


n 


0 




12022 


50 


0 


ooooo 


LM 


0 


0 






pO 


0 


00u;<j 


Itl 


) 


Cl 




13 


20 


1 


06500 




1 


TEMP 


CCCi’-iA /i.'i. • 




12 


0 


1215i- 






JJELT.'-. 






f5 


4 


600^ 0 


SlJ 


4 


aiWFL 






50 


0 


CX)000 


lil 


( , 
-i 


0 




1 'i 


-;o 


0 


•'LKOO 


. -r 

lu; X 


rt 

• j 


,1 






50 


0 


.-)OCV^T 




1 


0 
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12026 

12027 

12030 

12031 

12032 

12033 

1203h 

12035 

12036 

12037 

12040 



50 0 00000 

50 0 00000 
20 1 05500 
33 1 06500 

31 0 12175 REPL 

22 0 12141 
22 2 12141 

51 1 00001 
12 0 12152 
30 0 12151 
20 0 12151 
75 0 12017 

56 1 12105 NEXT 

56 1 12073 
56 1 12111 
56 1 12051 
56 1 12074 
56 1 12101 
56 1 12115 
50 1 00000 
12 1 06500 

32 1 06500 



ENI 


0 


0 


ENI 


0 


0 


STA 


1 


HOLD 


FDV 


1 


TEMP 


FSB 


0 


C»ICE 


AJP 


0 


FIRST 


MP 


2 


FIRST 


INI 


1 


1 


IDA 


0 


DELTA+1 


FAD 


0 


DELTA 


STA 


0 


DELTA 


SLJ 


0 


BOI 


SIU 


1 


JAY+3 


SIU 


1 


INCL+5 


SIU 


1 


JAYTO+3 


SIU 


1 


RMK 


SIU 


1 


INCL46 


SIU 


1 


JAY-1 


SIU 


1 


JAXFR+3 


ENI 


1 


0 


IDA 


1 


TEMP 


FMU 


1 


TEMP 



SINE TO STORAGE 
SINE/COSINE 
SINE/COSINE-R/H+T 
JUNPS IF A IS 0 
JUMPS IF A IS + 



STARTS FORMATION OF 
THE LECENDRE POL. 



COSB IN A REG. 
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I 



1 

I 



I 







* ^ 






1204l 

12042 

12043 

12044 

12045 

12046 

12047 

12050 

12051 

12052 

12053 

12054 



20 


0 


12173 




STA 


0 


DON 




32 


0 


12157 




FMU 


0 


FLT2 




31 


0 


12155 




FSB 


0 


FLT 




53 


0 


12156 




FDV 


0 


FLTl 


P2 12^ A REG. 


20 


1 


07500 




STA 


1 


W0RK2 




12 


0 


12173 




LDA 


0 


DON 




32 


0 


12161 




FMU 


0 


flt4 




20 


0 


12174 




STA 


0 


DON+1 




12 


0 


12173 




LDA 


0 


DON 




32 


0 


12173 




FMU 


0 


DOI 




32 


0 


12162 




FMU 


0 


FLT5 




31 


0 


12174 




FSB 


0 


DCN+1 




30 


0 


12157 




FAD 


0 


FLT2 


P4 IN THE A REG. 


33 


0 


12160 




FDV 


0 


FLT3 




20 


1 


10500 




STA 


1 


WOR}^ 




50 


0 


00000 




ENI 


0 


0 




54 


1 


00000 


RMK 


ISK 


1 


0 


»IEER OF AJ^GLES 


75 


0 


12040 




SU 


0 


NEXT+4 


INSERTED BY PROGRAM 


12 


0 


12164 


HALF 


LDA 


0 


THK 


FORMATION OF XB 


33 


1 


06500 




FDV 


1 


TEMP 




20 


1 


11500 




STA 


1 


HANS 




50 


0 


00000 




ENI 


0 


0 




54 


1 


00000 




ISK 


1 


0 


N1MER INSERTED 


75 


0 


12052 




SU 


0 


HALF 


m PROGRAM 
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I 

I 



f 

I 

I 

I 






fi 




1 



12055 

12056 

12057 

12060 

12061 

12062 

12063 

1206^1 

12065 

12066 

12067 



12 


1 


0 

0 

H 

H 




LDA 


1 


HMS 




32 


0 


12166 




FMJ 


0 


TAU 




20 


1 


11500 




STA 


1 


H/IJS 




50 


0 


00000 




ENI 


0 


0 




54 


1 


00000 




ISK 


1 


0 


NUMBER ENTERED 


75 


0 


12055 




SU 


0 


1^2 


W PROGRAM 


12 


1 


11500 


POWER 


LDA 


1 


HANS 


STARTS F0Rlt4ATI0N 


75 


4 


60600 




SU 


4 


EXPFL 


OF E TO THE TAU X 


50 


0 


00000 




ENI 


0 


0 




50 


0 


00000 




ENI 


0 


0 




20 


1 


12500 




STA 


1 


SOLD 




50 


0 


00000 




ENI 


0 


0 




54 


1 


00000 




ISK 


1 


0 




75 


0 


12060 




SU 


0 


POWER 




54 


2 


00001 




ISK 


2 


1 




75 


0 


12066 




SU 


0 


/+2 




75 


0 


12076 




SU 


0 


HANG 




50 


0 


00000 




ENI 


0 


0 




50 


1 


00000 


INCL 


ENI 


1 


0 




51 


1 


00001 




INI 


1 


1 




12 


0 


12165 




LDA 


0 


HGT 


STARTS CALCULATION 


33 


1 


06500 




FDV 


1 


TEMP 


R CSCB“cSECB 
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!; 



'I 

I 

I 



i 



12070 

12071 

12072 

12073 

1207^ 

12075 

12076 

12077 

12100 

12101 

12102 

12103 



20 


1 


12500 




STA 


1 


SOLD 






12 


0 


12163 




IDA 


0 


RAD 






33 


1 


05500 




FDV 


1 


HOLD 






31 


1 


12500 




FSB 


1 


SOLD 






32 


0 


12166 




FMU 


0 


TAU 






20 


1 


11500 




STA 


1 


HANS 






5** 


1 


00000 




ISK 


1 


0 






75 


0 


12067 




SLJ 


0 


/-4 






50 


1 


00000 




ENI 


1 


0 






56 


1 


12063 




SIU 


1 


POWER+3 






50 


1 


00000 




ENI 


1 


0 






75 


0 


12060 




SLJ 


0 


POWER 






12 


0 


12155 


HAKG 


IDA 


0 


FLT 






33 


1 


12500 




FDV 


1 


SOLD 






20 


1 


12500 




STA 


1 


SOLD 






12 


0 


12155 




LDA 


0 


FLT 






31 


1 


12500 




FSB 


1 


SOLD 






20 


1 


12500 




STA 


1 


SOLD 






$k 


1 


00000 




ISK 


1 


0 






75 


0 


12076 




SU 


0 


H/DG 






12 


0 


12153 


JPY 


IDA 


0 


WORK 


STARTS 


CALCULATION 


32 


1 


12500 






1 


SOLD 


OF THE 


DfTERGAL 


32 


1 


05500 




FMU 


1 


HOLD 






32 


0 


12152 




FMU 


0 


DELTA+1 
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I 



1 



12104 

12105 

1210b 

12107 

12110 

12111 

12112 

12113 

12114 

12115 

12116 



30 0 13500 
20 0 13500 
54 1 00000 
75 0 12102 
12 1 07500 JAYTO 
32 1 12500 
32 1 05500 
32 0 12152 
30 0 13504 
20 0 13504 
54 1 00000 
75 0 12106 
12 1 10500 JAYFR 
32 1 12500 
32 1 05500 

32 0 12152 
50 0 13510 
20 0 13510 
54 1 00000 
75 0 12112 
12 0 13504 

33 0 13500 



FAD 


0 


JAYO 


STA 


0 


JAYO 


ISK 


1 


0 


SU 


0 


/-3 


LDA 


1 


W0RK2 


FMJ 


1 


SOLD 


FMU 


1 


HOLD 


Fl^ 


0 


DELTA+1 


FAD 


0 


JAY2 


STA 


0 


JAY2 


ISK 


1 


0 


SU 


0 


/“3 


IDA 


1 


W0RK4 


FMU 


1 


SOLD 


FtlU 


1 


HOLD 


FMJ 


0 


DELTA+1 


FAD 


0 


JAY4 


STA 


0 


JAY4 


ISK 


1 


0 


SU 


0 


/“3 


LDA 


0 


JAY2 


FDV 


0 


JAYO 



KUNmER OF TIMES 
TIIX R/K 



R^J^ffiER OF TIMES 
TILL R/H 
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i 



12117 

12120 

12121 

12122 

12125 

1212U 

12125 

12126 

12127 

12150 

12131 

12152 



20 


O 

H 


STA 


0 


JAY20 


12 


0 15510 


LDA 


0 


JAY4 


55 


0 15500 


FDV 


0 


JAYO 


20 


0 15520 


STA 


0 


jay 4 o 


12 


0 1551^ 


LDA 


0 


JAY20 


52 


0 1551^ 


FMU 


0 


JAY20 


20 


0 15524 


STA 


0 


JAY4044 


52 


0 12167 


FMU 


0 


THEOR 


20 


0 15550 


STA 


0 


jay 4 o+io 


12 


0 15520 


LDA 


0 


jay 4 o 


52 


0 15520 


FMU 


0 


jay 4 o 


20 


0 15554 


STA 


0 


jay 4 o+i 4 


52 


0 12170 


FMJ 


0 


THEOR+1 


20 


0 15540 


STA 


0 


JAY4 o +20 


5^ 


1 77777 


ISK 


1 


77777 


75 


0 12126 


SLJ 


0 


/ 


50 


0 00000 


ENI 


0 


0 


7k 


7 42 o 4 o 


EXF 


7 


42 o 4 o 


75 


4 71000 


SU 


4 


71000 


50 


0 00000 


ENI 


0 


0 


01 


0 15500 


01 


0 


JAYO 


oh 


0 00001 


o4 


0 


1 


12 


0 12154 


LDA 


0 


MIKE 


70 


0 12151 


RAD 


0 


/-I 
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12133 

1213 *+ 

12135 

12136 

12137 

1214 o 

12li»l 

12142 

12143 

12144 

12145 



54 5 00010 


ISK 


5 


10 




75 0 12130 


SLJ 


0 


/-5 




54 1 77777 


ISK 


1 


77777 




75 0 12134 


SU 


0 


/ 




54 2 77777 


ISK 


2 


77777 




75 0 12135 


SU 


0 


/ 




50 0 00000 


ENI 


0 


0 




74 7 42o4o 


EXF 


7 


42 o 4 o 




74 0 42003 


EXF 


0 


42003 




50 0 00000 


ENI 


0 


0 




76 0 05000 


SLS 


0 


5000 




50 0 00000 


ENI 


0 


0 




12 0 12147 FIRST 


IDA 


0 


SEC 


MODIFIES THE PROGRAM 


20 0 12030 


STA 


0 


REPL 


WHEN THE VALUE OF 


12 0 12150 


IDA 


0 


SEC+1 


B CHAN®S 


20 0 12031 


STA 


0 


REPL+1 




56 1 12054 


SIU 


1 


HALF+2 




56 1 12057 


SIU 


1 


HAlF+5 




56 1 12063 


SIU 


1 


POWER+3 




56 1 12075 


SIU 


1 


INCL+7 




56 1 12066 


SIU 


1 


INCL 




51 1 00001 


INI 


1 


1 
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12146 

1214? 

12150 

12151 

12152 

12153 

12154 

12155 

12156 

12157 

12160 

1216 i 



75 


0 12017 




SLJ 


0 BGN 




50 


0 00000 




ENI 


0 0 




31 


0 12172 


SEC 


FSB 


0 PLEAS 


PARAMETER TO MODIFY 


22 


0 12034 




AJP 


0 NEXT 


PROGRAM FOR B 


22 


2 12034 




MP 


2 NEXT 


GRE.ATER THAN R/H+T 


51 


1 00001 




nil 


1 1 




00 


0 00000 


DELTA 


OCT 


0 




00 


0 00000 










17 


6 64 o 6 i 




DEC 


.001 




11 


5 64570 










20 


0 i4ooo 


WORiC 


DEC 


1.0 




00 


0 00000 










00 


0 00004 


MIKE 


0 


0 4 




00 


0 00001 




0 


0 1 




20 


0 i4ooo 


FLT 


DEC 


1.0 




00 


0 00000 










20 


0 24000 


FLTl 


DEC 


2.0 




00 


0 00000 










20 


0 26000 


FLT2 


DEC 


3.0 




00 


0 00000 










20 


0 44000 


FLT5 


DEC 


8.0 




00 


0 00000 










20 


0 57400 


flt4 


DEC 


30.0 





00 0 00000 
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12162 

12165 

12164 

12165 

12166 

12167 

121T0 

12171 

12172 

12175 

12175 



DEC 55.0 



20 0 64500 
00 0 00000 
20 0 25050 
75 5 ^1217 
20 0 25050 
75 3 *0217 
20 0 54400 
00 0 00000 
17 7 46215 
20 7 12601 
17 7 46416 
25 4 02050 
17 7 1*^511 
51 2 62577 
00 0 00000 
00 0 00000 
00 0 00000 
00 0 00000 
00 0 00000 
00 0 00000 
00 0 00000 
00 0 00000 



FLT 5 



RAD DEC 



THK DEC 



HGT DEC 



TAU DEC 



THEOR DEC 



DEC 



ZRO OCT 



PLEAS BSS 



DON BSS 



ONCE BSS 



2.5*t 

2.54 

18.0 

,098 

.1020 

.00907 

0 

1 

2 

1 



END 
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f 
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